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Abstract 

We have used the lowest order constrained variational (LOCV) method to calculate some ground 
state properties of polarized liquid s He at zero temperature with the spin-dependent correlation 
function employing the Lennard- Jones and Aziz pair potentials. We have seen that the total energy 
of polarized liquid s He increases by increasing polarization. For all polarizations, it is shown that 
the total energy in the spin-dependent case is lower than the spin- independent case. We have seen 
that the difference between the energies of spin-dependent and spin-independent cases decreases by 
increasing polarization. We have shown that the main contribution of the potential energy comes 
from the spin-triplet state. 
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INTRODUCTION 



Helium has two stable isotopes: that of mass 4 is readily available as helium gas or 
liquid from the atmosphere or gas wells, while that of mass 3 is extremely rare in nature 
and only became available commercially in the 1950s Liquid 3 He is particularly suited 
to study correlation among the strongly interacting many-body fermionic systems. Several 
approaches have been used for investigating the properties of normal liquid 3 He. These 



CBPT formalism 



a, 



are mainly based on the STLS scheme 2|, mott localization 3j, spin fluctuation theory 
Green's function Monte Carlo (GFMC) [sj , FN-DMC, DMC, VMC and EMC simulations 



nonperturbative renormalization group equation 



functional formalism jt| , correlated basis functions (CBF) [lOj and Fermi hyper-netted 



nonlocal density 



chain (FHNC) The spin polarized liquid 3 He as an interesting many-body system has 



been investigating different approaches snch as FHNC M , GFMC Q, CBF M and 
transport theory [15]. 

In recent years, we have studied both normal and polarized liquid 3 He at zero and finite 
temperature [if], 0, Q, Q] • In these calculations, the lowest order constrained variational 
(LOCV) method based on the cluster expansion of the energy functional has been used. 
This method is fully self-consistent, since it does not introduce any free parameter to the 
calculations. We have also used the LOCV method in many-body calculations of dense 



matter 



20j. Recently, we have used this method to calculate some properties of the polarized 



neutron matter and the polarized symmetrical and asymmetrical nuclear matters [21|]. In 
these works, a comparison of our results and those of other many-body techniques indicates 
that the LOCV method is a powerful microscopic technique to calculate the properties of 
the polarized matter. 

In this work, we use the LOCV method to compute the ground state energy of the polar- 
ized liquid 3 He at zero temperature b y e mploying; the spin-dependent correlation function 



with the Lennard- Jones 



22 1 and Aziz 



23 



24| pair potentials. 



2 



II. LOWEST ORDER CONSTRAINED VARIATIONAL METHOD 



We consider a system of N interacting 3 He atoms with N + spin up and iV spin down 
atoms. The total number density (p) and spin asymmetry parameter (£) are defined as 



P 



p + + p , 

N + - N~ 
N 



(1) 



£ shows the spin ordering of matter which can have a value in the range of £ = 0.0 (unpo- 
larized matter) to £ = 1.0 (fully polarized matter). For this system, we consider the energy 
per particle up to the two-body term in the cluster expansion, 



E = E 1 + E : 



2- 



(2) 



where 



Ex = ^(37^)31(1 + 0* + (1-0*1, 



(3) 



In the above equation, w(12) is the effective pair potential, 

w(12) = ~[F{12), [V? 2 , F(12)]] + F(12)V(12)F(12), 



(4) 



where E(12) is the two-body correlation operator and V(12) is the pair potential between 



the helium atoms. In our calculations, we use the Lennard- Jones 
potentials. The Lennard- Jones pair potential is as follows 



2j, 



22] and Aziz 



V(r) = 4e 



12 



24( | pair 



(5) 



where 



e = 10.22ZT, a = 2.556A • 



The Aziz pair potential has the following form [23| 
V(r) = e^Ae~ ar/rm - 

where 



r \ 6 / r \ 8 / r \ 10 

c 6 (-) + c 8 p) + c 10 p 

r / V r / V r 



F(r) 



(6) 



(7) 



F(r) 



1 ; — > D, 

I TV™ ' 



(8) 



and 



= 10.8K, A = 0.5448504 x 10 b , 
a = 13.353384, r m = 2.9673A, 
C 6 = 1.37732412, C s = 0.4253785, 



(9) 



Cio = 0.178100, D = 1.241314- 
Now, we consider a spin-dependent correlation function as follows 

F(12) = /o(r 12 )P + /i(ri 2 )Pi, (10) 

where 

Po = ^(1 - ffi.<T 2 ), 

Pi = ^(3 + (7i.(7 a ). (11) 

/o and /i indicate the spin-singlet and spin-triplet two-body correlation functions, respec- 
tively. With the above two-body correlation function, we have derived the following relation 
for the effective pair potential, 

w s {r) = ^-{fJ{r)f + f s {r)V{r), (12) 

and then the two-body energy E 2 is found by 

POO 

E 2 = 2np V" / drr 2 w s (r)a s . (13) 
s=0,l Jo 



InEq. ([T3D 

*i = \[(l + e(l-l 2 (k F+ r)) + (l-0 2 (l-l 2 (k F -r)) 

+(1 l(k F+ r)l(k F -r))}. (14) 

k F ± = (Qir 2 p ± )^ is the Fermi momentum and l{x) is given by 

3 

l(x) = — [sin(x) — xcos(x)]. (15) 
Now, we minimize the two-body energy Eq. f|T3|) with respect to the variations in the 



two-body correlation function subject to the normalization constraint 

1 



25] 



2>1tf(12) -F 2 (12)|y-ji> = l- (16) 



The normalization constraint is conveniently re-written in the integral form as 

4tt P £ / drr 2 [h 2 (r) - f*(r))a s = 1, (17) 

where the Pauli function h(r) is 

h(r) = {1 - + 2 l 2 (k F+ r) + (1 - 2 / 2 (^-r)]}-i (18) 

The minimization of the two body energy E 2 gives the following Euler-Lagrange differential 
equation for the two-body correlation function f s (r), 

/;'(r) + (- + %» - ™ (V(r) - 2A) = 0. (19) 
r a s n 

The Lagrange multiplier A imposes by normalization constraint. For s = and s = 1 states, 
the two-body correlation function / s (r) is obtained by numerically integrating Eq. (jl~9l) . 
Using this two-body correlation function we can determine the effective pair potential w s (r) 
as a function of interatomic distance from Eq. (|T2|) . Finally, the two-body energy E 2 and 
the total energy of system can be calculated. 



III. RESULTS AND DISCUSSION 



We have calculated some ground state properties o 



23, 



' the polarized liquid 3 He at zero 



24| pair potentials using the spin- 



temperature with the Lennard- Jones [22| and Aziz 
dependent correlation function. Our results £11*6 cLS follows. 

The spin-dependent two-body correlation functions at s = state and s = 1 states for 
different values of spin asymmetry parameter (£) are shown in Fig. [Q These figures show 
that the correlation function at s = 1 state heals to pauli function, h(r), more rapidly than 
s = state. Therefore, the s = 1 state has a shorter correlation length with respect to 
s = state. For large values of r, fo(r) and fi{r) have the same values and therefore, the 
spin-dependent part of correlation operator (Eq. [TO]) is vanished. In these figures, the spin- 
independent two-body correlation function are also plotted for comparison. It is seen that 
the spin-dependent correlation function differs from spin-independent correlation function, 
except for the fully polarized matter (£ = 1.0). This is due to the fact that for fully polarized 
matter, there is only s = 1 state. From Fig. [lj we can see that the correlation functions 
with the Lennard- Jones and Aziz potentials are nearly identical. 



In Fig. [21 we have shown the total energy of polarized liquid 3 He versus number density 
calculated both with the spin-dependent correlation and the spin-independent correlation at 
different values of spin asymmetry parameter £. We can see that the total energy increases 
by increasing £. Fig. [2] indicates that in the spin-dependent case, the total energy of the 
liquid 3 He is lower than the spin- independent case. It is also seen that for all values of £, the 
energy curve has a minimum which shows the existence of a bound state for this system. It is 
shown that the difference between the energies of spin-dependent case and spin-independent 
case decreases by increasing £ and it becomes zero as £ approaches to one. It is seen that 
for all values of the density and spin asymmetry parameter, the total energy with the Aziz 
pair potential is grater than that of the Lennard- Jones pair potential. 

The potential energy of the polarized liquid 3 He for different values of £ are presented in 
Fig. [31 for spin-dependent and spin-independent cases. This figure indicates that the poten- 
tial energy decreases by increasing the polarization. According to the above results, we can 
conclude that the increasing of kinetic energy dominates and this leads to the increasing of 
total energy by increasing £. Fig. [3] shows that the potential energy in the spin-dependent 
case has lower values with respect to the spin-independent case. It is also seen that the differ- 
ence between the potential energies of spin-dependent and spin-independent cases decrease 
by increasing £. We see that the potential energies with the Aziz and the Lennard- Jones 
pair potentials are different. This difference increases by increasing the density. 

In Fig HI the potential energies of s = and s — 1 states for different values of £ are 
compared. We have seen that the potential energy at s = 1 state is lower than at s = state. 
It can be concluded that the spin-triplet state has the main contribution in the potential 
energy of polarized liquid 3 He. It is also seen that the potential energy of s = (s = 1) 
state increases (decreases) by increasing £. For s = state, we can see that at low densities, 
the potential energies with the Lennard- Jones and Aziz pair potentials are nearly identical. 
However, for s = 1 state and high densities, the difference between these potential energies 
becomes appreciable. 

The equation of state of polarized liquid 3 He, P(p,£), can be obtained using 

Pfcfl^^M (20) 

In Fig. El we have presented the pressure of liquid 3 He as a function of the density (p) for 
fully polarized (£ = 1.0) and unpolarized (£ = 0.0) cases. This figure shows that for different 
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values of the polarization, the equations of state of liquid 3 He are nearly identical. From 
Fig. [5l it is seen that for both £ = 1.0 and £ = 0.0, the equation of state with the Aziz pair 
potential is stiffer than that of the Lennard- Jones pair potential. 

IV. SUMMARY AND CONCLUSION 

We have considered a system consisting of Helium atoms ( 3 He) with asymmetrical spin 
configuration and derived the two-body term in the cluster expansion of the energy func- 
tional by employing spin-dependent correlation function. Then, we have minimized the 
two-body energy term under the normalization constraint and obtained the Euler-Lagrange 
differential equation. By numerically solving this differential equation, we have computed 
the correlation function and then calculated the other properties of this system with the 
Lennard- Jones and Aziz pair potentials. It is shown that for the two different spin-singlet 
and spin-triplet states, the correlation functions are different from each other. Our results 
show that the introduction of the spin-dependent term in the correlation operator reduces 
the total energy of system by about 10%. It is also shown that the total energy increases 
by increasing the polarization. The difference between the energies of the spin-dependent 
and spin-independent cases decreases by increasing the polarization. We have seen that, 
the potential energy of these states have a remarkable difference. It is shown that the main 
contribution of the potential energy comes from s — 1 state. Our calculations show that 
there is a difference between the results with the Lennard- Jones and Aziz pair potentials, 
especially at high densities. 
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r(A) r(A) 



FIG. 1: The correlation function with the Aziz (dashed curves) and Lennard- Jones (full curves) 
pair potentials in the case of spin-dependent at s = and s = 1 states for £ = 0.0 (a), £ = 0.33 
(b), £ = 0.66 (c) and £ = 1.0 (d). Our results for the spin-independent correlation function are 
also presented for comparison. 
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Density (A 3 ) Density (A 3 ) 

FIG. 2: Our results for the total energy of the polarized liquid 3 He with the Aziz and Lennard- 
Jones (LJ) pair potentials in the case of spin-dependent (full curve) and spin-independent (dotted 
curve) correlation functions for £ = 0.0 (a), £ = 0.33 (b), £ = 0.66 (c) and £ = 1.0 (d). 
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FIG. 4: Our results for the potential energy of the polarized liquid 3 He with the Aziz and Lennard- 
Jones (LJ) pair potentials at s = (dotted curve) and s = 1 (full curve) states for £ = 0.0 (a), 
£ = 0.33 (b), f = 0.66 (c) and f = 1.0 (d). 
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FIG. 5: The equation of state of the fully polarized (dotted curve) and unpolarized (full curve) 
liquid s He with the Aziz and Lennard-Jones (LJ) pair potentials. 
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